Heat conduction in one dimensional nonintegrable systems 
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Two classes of ID nonintegrable systems represented by the Fermi- Pasta-Ulam (FPU) model and 
the discrete 4> 4 model are studied to seek a generic mechanism of energy transport in microscopic level 
sustaining macroscopic behaviors. The results enable us to understand why the class represented 
by the 4> 4 model has a normal thermal conductivity and the class represented by the FPU model 
does not even though the temperature gradient can be established. 
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Heat conduction in one-dimensional (ID) noninte- 
grable Hamiltonian systems is a vivid example for study- 
ing microscopic origin of the macroscopic irreversibility 
in terms of deterministic chaos. It is one of the oldest but 
a rather fundamental problem in nonequilibrium statisti- 
cal mechanics Q . Intended to understand the underlying 
mechanism of the Fourier heat conduction law, the study 
of heat conduction has attracted increasing attention in 
recent years 

Based on previous studies, we can classify the ID lat- 
tices into three categories. The first one consists of inte- 
grable systems such as the harmonic chain. It was rig- 
orously shown jl2j that, in this category no tempera- 
ture gradient can be formed, and the thermal conductiv- 
ity is divergent. The second category includes a num- 
ber of nonintegrable systems such as the Lorentz gas 
model [p|JlC|l, the ding-a-ling and alike models and 
the Frenkcl-Kontorova (FK) model |J and etc.. In this 
category, the heat current is proportional to iV -1 and 
the temperature gradient dT/dx ~ N~ 4 , thus the ther- 
mal conductivity n is a constant independent of system 
size N. The Fourier heat conduction law (J = —ndT/dx) 
is justified. The third category also includes some non- 
integrable systems such as the FPU |||| chain, the di- 
atomic Toda chain 0] , the (mass) disorder chain || , and 
the Heisenberg spin chain |)| and so on. In this cate- 
gory, although the temperature gradient can be set up 
with dT/dx ~ N~ 4 , the heat current is proportional to 
N a ~ 1 with a ~ 0.43, the thermal conductivity k ~ N a 
which is divergent as one goes to the thermodynamic 
limit N — > oo. 

These facts suggest that the nonintegrability is neces- 
sary to have a temperature gradient, but it is not suffi- 
cient to guarantee the normal thermal conductivity in a 
ID lattice. This picture brings us to ask two questions 



of fundamental importance: (i) Why do some noninte- 
grable systems have normal thermal conductivity, while 
the others fail? (ii) How can the temperature gradient 
be established in those nonintegrable systems having di- 
vergent thermal conductivity? 

The reason for the divergent thermal conductivity in 
integrable system is that the energy transports freely 
along the chain without any loss so that no temperature 
gradient can be established. The set up of temperature 
gradient in nonintegrable systems implies the existence 
of scattering. However, the different heat conduction be- 
havior in the two categories of nonintegrable systems in- 
dicates that the underlying mechanism must be different. 
To get the point, let's write the Hamiltonian of a generic 
ID lattice as 



Hi = 



EL 

2 



+ V{x l - 1 ,x i ) + U(x t ), (1) 



where V(xi-i, Xi) stands for the interaction potential of 
the nearest-neighbor particle and U{xi) is an external 
(or on-site) potential. The origin of external potential in 
real physical systems varies from model to model. For 
instance, in the FK model || the external potential is 
the interaction of the adsorbed atoms with the crystal 
surface. It is U(x) that distinguishs from the two cate- 
gories of nonintegrable lattices. U (x) vanishes in all ID 
lattices having divergent thermal conductivity. We are 
thus convinced to conclude that the external potential 
plays a determinant role for normal thermal conduction. 

In this paper we would like to study the scattering 
mechanism and the role of the external potential in 
heat conduction in the two categories of nonintegrable 
systems. For this purpose, we choose two representa- 
tives from these two categories, i.e. the discrete cj) 4 
model (see, e.g. Ref. ]U|) and the FPU model. Both 
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models are the simplest anharmonic approximation of a 
monoatomic solid. In the 4> 4 model, V takes the harmonic 
form, and the external potential U(x) = mx 2 /2 + fix 4 /A, 
with m fixed to be zero in this paper. In the FPU 
model, U vanishes and V takes the anharmonic form of 
(xi — Xi-i) 2 /2 + (3(xi — £i_i) 4 /4, and (3 — 1 throughout 
this paper. In the case of fi = 0, the FPU model reduces 
to the harmonic chain. 

In our numerical simulations the Nose-Hoover ther- 
mostates jlij are put on the first and the last particles, 
keeping them at temperature T + and T_, respectively. 
The motion of these two particles are governed by 

Xi = ~(+Xi + fi - / 2 , C+ = il/T+ 



1 



x N 



<-x N + f N -f N+1 , C- = ii/T--l. (2) 



where fi = —(V + U') is the force acting on the i'th 
particle. The equation of motion of other particle is 
Xi = fi — fi+i. The eighth-order Runge-Kutta algorithm 
was used. All computations are carried out in double 
precision. Usually the stationary state set in after 10 7 
time units. We should point out that we have performed 
computations by using other types of thermostate, and 
no qualitative difference has been found. 
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(a) Temperature profiles for the FPU model, the 





FIG. 1 

4> 4 model, the harmonic and the monoatomic Toda model, 
(b) The quantity J x N versus the system size N for the 
FPU (solid square) and the 4 (solid circle) . That of the har- 
monic and the monoatomic Toda are shown in the inset of 
(b). The lines in (b) and its inset are drawn for guiding the 
eye. T+ = 0.3, T_ = 0.2. 

Fig. 1(a) shows temperature profiles. In all noninte- 
grable systems, the temperature scales as T — T(i/N). 
However, in the FPU case there is a singular behavior 
near the two ends, which is a typical character of ID 
nonlinear lattices having divergent thermal conductivity. 
In the same figure we also show the temperature pro- 
files for two integrable lattices: the harmonic and the 
monoatomic Toda models. In these two cases no temper- 
ature gradient could be set up and the stationary state 
corresponds to T = (T+ +T_ )/2, which is consistent with 
the rigorous result pjfl. 



In Fig. 1(b), we plot the quantity J x N versus 
N for the FPU model and the </> 4 model. The inset 
shows the same quantity for the harmonic chain and the 
monoatomic Toda chain. The local heat flux is defined 
by Ji — it a ® V +1 ■ We found that when the system reaches 
a stationary state, the time average (Jj) is site indepen- 
dent, it is denoted as J. For the harmonic chain and the 
monoatomic Toda chain J x N is expected to be propor- 
tional to N since J is TV-independent. This is indeed the 
case as illustrated in the inset. In both the FPU and the 
4> A models dT/dx is proportional to —1JN, the thermal 
conductivity k = -J/{dT/dx) cx J x N. Fig. 1(b) tells 
us that, in contrast to the integrable systems and the 
FPU model, heat conduction in the </> 4 model obeys the 
Fourier law. 
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FIG. 2. The momentum excitation in different lattices, 
(a) the harmonic chain; (b) the monoatomic Toda chain; (c) 
the FPU model; and (d) the </> 4 model. Inset in (d) illustrates 
the time evolution of the total momentum of the d> 4 chain. 



The heat current J in all nonintegrable systems de- 
creases as the system size N is increased (J ~ iV" -1 , < 
a < 1). To clarify the underlying mechanism we decom- 
pose the interaction of the thermostat into a series of 
kicks and study the transport of a single kick along the 
chain. A free boundary condition is used in our calcula- 
tion, but we should stress that the results do not depend 
on the types of boundary condition. In Fig. 2 we plot pi 
versus i after a long time (t = 800) for four lattices: the 
harmonic chain (a); the monoatomic Toda chain (b); the 
FPU model (c); and the </> 4 model (d). The amplitude of 
the wave profile in the harmonic chain decreases continu- 
ously with time, but the global profile keeps unchanged. 
In both the FPU and the Toda models, we observe a 
solitary wave separates from the long tail. Initially this 
wave front is connected with other low amplitude excita- 
tions. After a certain time, this wave front, moves faster, 
separates from the tail and goes forward and keeps its 
amplitude. In the mean time the tails behind it evolve 
in the same way as that in the harmonic chain. In the 
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</> 4 model, the head part of the profile becomes weaker 
and weaker. The reason is that in the first three cases 
(a-c) both the total energy and the total momentum are 
conserved, whereas in the cf> 4 model the momentum con- 
servation breaks down due to the external potential. The 
inset in Fig. 2(d) shows that the total momentum in 
the </> 4 model decreases at least exponentially with time. 
The decay of the momentum with time indicates a loss 
of correlation. It is thus reasonable to envisage the en- 
ergy transport along the cf) A chain as a random walk-like 
scattering. 

The solitary waves in the FPU chain exchange energy 
and momentum when colliding with each other. It causes 
the energy loss, and the heat current decreases when the 
system size is increased. To show this, we start two ex- 
citations at the two ends of the chain with different mo- 
mentum, one moves to the right and another left. Let 
Pi = 6, pn = 3 and = 0, i ^ 1 and N as our initial 
excitations. We calculate the momentums of the solitary 
waves (by simply summing up momentums of several lat- 
tices around the peaks) and investigate its change before 
and after the interaction. We find that the bigger one 
generally transfers part of its momentum and energy to 
the smaller one, as is shown in Fig. 3(a). The collision 
takes place at t = 850, where a peak is shown up. 
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FIG. 3. (a) The time dependence of the momentum of a 
solitary wave, (b) The momentum (after collision) versus the 
phase. Solid triangle represents the results of the case in which 
both left and right solitary waves have the same initial mo- 
mentum. Solid circle represents the case of the solitary waves 
having different initial momentum (see text for more in de- 
tail). The horizontal line is the momentum before collision, 
(c) The maximal momentum gain Ap max versus initial mo- 
mentum p . (d) The ratio of the heat flux of J(400)/ J(800) 
versus the average temperature T — (T + + T_)/2 in the 
semi-logarithmic scale. The results shown here are for the 
FPU model. 

Moreover, the interaction between solitary waves is 
found to depend closely on a "phase" difference. Here 
the "phase" difference is defined as a time lag between 



the excitations of two solitary waves. For instance, if we 
excite a solitary wave from the left end at time t, and 
another one from the right end at time t + 6, then 5 is 
the "phase" difference. These two solitary waves, trav- 
eling through the chain in opposite directions, will col- 
lide with each other after a certain time. Although the 
physical meaning of the "phase" is not obvious, it is an 
important and good quantity to describe the interaction. 
We show pi versus S for two different kinds of collision 
in Fig 3(b), where pi is the momentum of the solitary 
wave from the left after collision. In the first case both 
left and right solitary waves have the same initial mo- 
mentum p^ — pji = 6.27, which is excited by an initial 
condition of p\ = Pn = 6 and pi = for other i. In 
the second case, the left one has pl = 6.27 and the right 
one pn — 3.28, excited by an initial condition of p\ = 6, 
Pn = 3 and pi = for i ^ 1 and N. The figure shows 
that P£ depends on the "phase" sinusoidally. 

Other interesting features of the collision of solitary 
waves are shown in Fig. 3(c), where we plot the maxi- 
mum momentum gain Ap max versus the initial momen- 
tum po for the FPU model. Ap max is measured by sub- 
tracting the initial momentum po from the maximum pl 
in Fig. 3b. First of all, this picture tells us that the ex- 
change of momentum and energy depends on the initial 
momentum and energy. Secondly, there exists a criti- 
cal momentum below which no energy exchange can take 
place. The critical ~ 1.8 is clearly seen in the figure. 
Forpo < Po, Ap max is zero. This result is very significant, 
it indicates that there exists a threshold for the solitary 
wave interaction, below this threshold the interaction 
ceases, i.e. no momentum and energy exchanges between 
the solitary waves. A direct consequence of this fact is 
the existence of a threshold temperature below which the 
FPU chain should behave like a harmonic chain, namely, 
the excited waves travel freely along the chain without 
any energy loss, no temperature gradient can be set up, 
and the heat current remains a constant even though the 
size of the chain is changed. To testify this argument, we 
show the quantity J(400) / J(800) versus T = (T + +T_)/2 
in Fig. 3(d), where J(N) is the heat current flux for 
the system of size N. In the case of a size-independent 
J(N) one should get J(400)/J(800) = 1, otherwise one 
would get J(400)/J(800) > 1. Fig. 3(d) captures this 
transition nicely for the FPU chain. The corresponding 
temperature threshold is about T c ~ 0.01. In the region 
of T ~ 0.001 the numerical calculations do show that no 
temperature gradient is formed. 

The different scattering mechanism in the FPU chain 
and those chains having normal thermal conductivity 
leads to a different temperature dependence of ther- 
mal conductivity k(T). In Fig. 4 we plot k(T) for a 
FPU chain with an external potential of form U(x) = 
— "fCos{2irx) for four different values of 7 = 0,0.01,0.05 
and 0.1. The chain size is fixed at N = 100. As pointed 
out above, in small 7 regime such as 7 = and 0.01, the 
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energy transport is assisted by the solitary waves, the 
system has a large k which decreases as the temperature 
is increased. However, in the opposite regime (7 = 0.05 
and 0.1), the energy transport is diffusive and obeys the 
Forier law, n increases with temperature, because more 
phonons are excited. 
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FIG. 4. The thermal conductivity «(T) for the FPU chain 
with an external potential U(x) — cos(2ttx) . 

In summary, by studying two classes of ID nonintc- 
grable lattices, we have answered the two questions raised 
at the beginning: (i) The multiple scattering of the ex- 
cited modes by the external potential leads to a decay 
(at least exponentially) of correlation, so that a diffusive 
transport process can be reached, and the heat conduc- 
tion obeys the Fourier law; (ii) Although the interaction 
of solitary waves makes it possible to set up temperature 
gradient in the FPU and alike nonintegrable models, the 
momentum conservation prohibits the diffusive transport 
and consequently leading to the divergent thermal con- 
ductivity. In addition, we have uncovered an important 
fact in the FPU model, namely, the existence of a thresh- 
old temperature, below which the FPU mode behaves like 
a harmonic chain. 
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sions. This work was supported in part by Hong Kong 
Research Grant Council and the Hong Kong Baptist Uni- 
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Note added in proof. - After submission of this pa- 
per we get to know the following results. Prosen and 
Campbell fl5|| proved in a more rigorous way that for 
a ID classical many-body lattice total momentum con- 
servation implies anomalous conductivity. The normal 
thermal conductivity in the 4 lattice has also been ob- 
served by Aoki and Kusnezov |ll| . The role of the exter- 
nal potential has been further studied by Tsironis et al 
E7|. 
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